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Abstract 

Stellar rotation depends on different parameters such as age, mass, initial chemical 
composition, initial angular momentum, and environment characteristics. The range 
of values of these parameters causes the dispersion in the rotation period distribu¬ 
tions observed in young stellar clusters/associations. We focus our investigation on 
the effects of different circumstellar environments on stellar rotation. More specifi¬ 
cally, we consider the effects of a perturber stellar companion on the accretion-disc 
lifetime at early evolution stages. 

We are searching in stellar Associations for visual triple systems where all stellar 
parameters are similar, with the only exceptions of the unknown initial rotation 
period, and of the circum-stellar environment, in the sense that one of the two 
about equal-mass components has a close-by third ’perturber’ component. 

In the present study we analyse the 35-Myr old visual triple system TYC 9300- 
0891-lAB -|- TYC 9300-0529-1 in the young Octans stellar association consisting of 
three equal-mass KOV components. We collected from the literature all information 
that allowed us to infer that the three components are actually physically bound 
forming a triple system and are members of the Octans Association. We collected 
broad-band photometric timeseries in two observation seasons. We discovered that 
all the components are variable, magnetically active, and from periodogram analysis 
we found the unresolved components TYC 9300-0891-1AB to have a rotation period 
P = 1.383 d and TYC 9300-0529-1 a rotation period P = 1.634 d. 

TYC 9300-0891-lA, TYC 9300-0891-lB, and TYC 9300-0529-1 have same masses, 
ages, and initial chemical compositions. The relatively small 16% rotation period 
difference measured by us indicates that all components had similar initial rotation 
periods and disc lifetimes, and the separation of 157 AU between the component 
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A and the ’perturber’ component B (or vice-versa) has been sufficiently large to 
prevent any significant perturbation/shortening of the accretion-disc lifetime. 


Key words: Stars: activity - Stars: low-mass - Stars: rotation - Stars: starspots - 
Stars: pre main sequence: individual : TYC 9300-0891-1, TYC 9300-0529-1 


1 Introduction 


The early stage of star’s life is characterized by the presence of an accretion 
disc that is magnetically coupled with the central star (see, e.g., Menard & 


Bertout 1999). This coupling implies complex exchanges of angular momen¬ 
tum between the disc and the central star (see, e.g., Bouvier et al. 2007). 
On one hand, the star gains angular momentum from the disc, on the other 
hand this excess angular momentum must be somehow dissipated. In fact, 
observations tell us that, until the star-disc interaction is effective, the star’s 


rotation period in most cases remains about constant 

Bouvier et al. 

1993 

Edwards et al.||1993 

Rebull et al.||2004 

Ingleby et al.||201 

4). This is commonly 

referred to as star-disc locking. Accretion-driven winds (A 

datt & Pudritz||2005 


Matt & Pudritz 

2008a 

Matt & Pudritz 

2008b 

) and mass ejection episodes 

caused by magnetospheric reconnection i 

Zanni & Ferreira 

2013 

) are among 


the possible mechanisms that have been proposed to remove the excess an¬ 
gular momentum. Meanwhile, other complex mechanisms, such as hydromag- 
netic instabilities and intense magnetic fields, also contribute to redistribute 
internal angular momentum to the star’s external envelope (see, e.g., Spada 


et al. 2010, Spada et al. 2011). After the inner disk dispersion, the stellar ra¬ 


dius continues to rapidly decrease while the star is approaching the Zero Age 
Main Sequence and the rotation rate spins up, despite the angular momen¬ 


tum losses due to magnetized stellar winds (Matt & Pudritz 2007). The disc 


lifetime is variable, of the order of a few Myr, but rarely longer than about 10 


Myr (Ribas et al. 2014 Ingleby et al. 2014). A variable disc lifetime implies 


variable duration of the star-disc locking phase. Indeed, the rotation rate of a 
PMS star at a given age depends on the accretion disc lifetime, in the sense 
that the shorter the disc lifetime the faster the rotation rate with respect to 
a counterpart star with longer disc lifetime. Therefore, the rotation rate can 
be used to derive information on the disc lifetime. However, the rotation also 
depends on the initial rotation period, on the mass and disc mass, and on 
the environment properties. We intend to investigate how the environment, 
and specifically the presence of a nearby companion, can shorten the disc life- 
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time and, consequently, make the rotation spin up to start earlier. Multiple 
stellar systems where the components a have equal age, initial chemical com¬ 
position, and equal mass are best suited for this purpose. In these systems, 
differences between the rotation periods of the components are expected to 
depend only on the initial rotation periods and disc lifetimes. Specihcally, we 
investigate the hypothesis according to which a close companion can enhance 
the disc gravitational dispersal allowing one component to start spinning up 
earlier that the other components. In our hrst of such cases, the triple system 
BD—211074 in the (3 Pictoris Association ( Messina et ah 1 2014), we found that 
to explain the observed difference of rotation periods between the primary A 
(Prot = 9.3 d) and the secondary B component (Prot = 5.4 d), a signihcant 
difference of lifetimes of their discs has to be invoked. Specihcally, the lifetime 
of the disc of component B has undergone a shortening owing to enhanced 
dispersal by dynamical effects from the close-by component C at 15 AU. Such 
disc dispersal has shortened the disc-locking phase, allowing component B to 
start its spinning up earlier than component A. 

However, the most promising approach to distinguish between different initial 
rotation rates and different disc locking timescales on a hrm basis is a statis¬ 
tical study based on a large number of cases. We are building such a sample 
by identifying suited stellar systems and deriving the rotation periods of their 
components. 

In the present paper, we present the results for another of such systems that 
has been studied as part of the RACE-OC project (Rotation and Activity 


Evolution in Open Clusters; Messina 2007 Messina et ah 2008, Messina et 


al. 2010), while measuring the rotation periods of the members of the about 


30-40 Myr old Octans association (Messina et al. 2011). 

This is the visual triple system TYC 9300-0891-lAB -|- TYC 9300-0529-1. 
This system consists of three equal-masses (0.95 Mq), ages (~35 Myr), and 
chemical compositions components physically bound that differ only for their 
circum-stellar environments and, in principle, for their initial angular mo¬ 
menta. Whereas TYC 9300-0529-1 is quite distant (~3600AU) from TYC 
9300-0891-lAB, on the contrary TYC 9300-0891-lA has a close-by (~157 AU) 
companion TYC 9300-0891-lB that may have acted as perturber of its (by 
now dispersed) accretion disc. Any difference we may hnd in their rotation pe¬ 
riods, should be ascribed to either the effect of the perturber on the accretion 
disc, and on its lifetime, or on different initial rotation periods. 

To measure the rotation periods of TYC 9300-0891-1 and TYC 9300-0529-1 
we carried out a dedicated photometric monitoring in two different observation 
seasons. In Sect. 2, we present the available information from the literature 
on these targets. In Sect. 3, we describe our photometric observations. Their 
analysis is presented in Sect.4, and 5, discussion of the results and conclusions 
are given in Sect. 6 and 7. 
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2 The triple system 


TYC 9300-0891-lAB and TYC 9300-0529-1 form a visual triple system, at a 
mean distance of 174 pc from the Sun, consisting of three physically bound 
components having equal mass (~ 0.95 Mq), KOV spectral type, and an age 
of about 35 Myr. 


TYC 9300-0891-lAB was spatially resolved into a close visual binary by Tycho 
(ESA 1997) that measured an angular separation p = 0.9", a position angle 
PA = 128.4° at the epoch 1991.25, and an integrated Johnson V magnitude 

The A and B components are 


Vab = 10.95±0.07 mag (Hog et ah 2000) 


listed as a physical pair with identical proper motions in the Tycho catalogue. 
As shown in Sect. 3, we measured a similar angular separation p = 0.95" and 
similar position angle PA = 135±10°. 

TYC 9300-0529-1 is at an angular distance p = 20.5" from TYC 9300-0891- 
lAB (ESA 1997) and its Johnson V magnitude (|Hog et ah 
11.75±0.11 mag. 


2000) is V = 


In the WDS catalog (The Washington Visual Double Star Catalog, Mason et 


al. 2001), TYC 9300-0891-lAB and TYC 9300-0529-1 are reported to be a 


physical pair on the basis of their common proper motions. We have collected 
from the literature the available kinematic information. Specihcally, we re¬ 


trieved from Tycho (ESA 1997), WDS (Mason et al. 2001), UCAC2 (Zacharias 


et al. 2004), and PPMXL ( Roeeser et al.[ 2010) catalogues the proper motion 
measurements, and from Torres et al. (2006) and Elliott et al. (2014) the ra¬ 


dial velocity (RV) measurements. Both proper motions and RV measurements 
of TYC 9300-0891-lAB and TYC 9300-0529-1 are undistinguishable within 
the uncertainties (see Table|^, indicating that TYC 9300-0891-lAB and TYC 
9300-0529-1 are physically bound and forms a triple system. 


Torres et al. (2008) found this system to be member of the young Octans 


stellar association. However, for no member of this association trigonometric 
parallaxes exist and kinematic distances are inferred. Recently, |Murphy~^ 


Lawson (2014) identihed 29 new low-mass members of Octans. Their study 


allowed to compute accurate mean space motion, the proper motions, radial 
velocities and, for the hrst time, a reliable age estimate of 30-40 Myr based 
on the Lithium Depletion Boundary. In their study, they inferred kinematic 
distances of 173 and 175 pc for TYC 9300-0891-lAB and TYC 9300-0529-1. 
We note that the difference between the radial velocities of TYC 9300-0891- 
lAB and TYC 9300-0529-1 is a factor 2.5 smaller that the radial velocity 


^ In the Tycho Double Star Catalogue (TDSC) ( Fabriclus et al.j2002 ) the following 
magnitudes are reported for both components Va = 11.50±0.07 mag and Vb = 
11.83±0.11 mag. 
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dispersion measured among the Octans members by Murphy & Lawson (2014), 


as well the difference between the proper motions is a factor 2 smaller that the 


proper motion dispersion measured among the Octans (see again Murphy & 


Lawson 2014). These circumstances further support that the three components 


are physical bound. 


The primary A component is a KOVe emission-line star (Torres et ah 2006), 


whereas no spectral classification was provided for the secondary B compo¬ 
nent. However, considering that they have similar V magnitudes (see Sect. 3), 
we can infer for the B component the same KOV spectral type. TYC 9300- 
0529-1 also has KOV spectral type (Torres et al. 2006). Therefore, we deal 


with a system with three about identical components. We note that [Torres et 
al. (2006) collected only one spectrum of TYC 9300-0891-lAB (unresolved), 
which exhibited emission lines, and one of TYC 9300-0529-1. However, we 
know that in the case of young stars, the emission features have variable EW 
and their detection is rotation-phase dependent (see, e.g. Parihar et al. 2009). 


Stars that show emission lines at some epochs may not at other epochs, de¬ 
pending on the activity level and the active region visibility along the line 
of sight at the epochs of observation. Additional spectra may reveal emission 
lines also in TYC 9300-0529-1. 


The system was detected by ROSAT as X-ray source (IRXS J184950.1-715703; 


Voges et al. 1999), although the components could not be resolved. Torres et al. 
(2006) measured Li EW = 310mA for the unresolved component TYC 9300- 


0891-lAB and 300 mA for TYC 9300-0529-1. From these values, da Silva et 


al. (2009) derived the Li abundances A(Li) = 2.69 and A(Li) = 2.64, adopting 
Teff = 5168 K and T^s = 5140 K for TYC 9300-0891-lAB and TYC 9300- 
0529-1, respectively. 


Torres et al. (2006) also measured the projected rotational velocity of TYC 
9300-0891-lAB {v sini = 9 kms“^) and TYC 9300-0529-1 {v sinz = 25.8 kms” 




Since all three components have late spectral type and fast rotation (see 
Sect. 4), all are expected to to be spotted and to exhibit photometric variability 
induced by such surface inhomogeneities. Actually, our own magnitude mea¬ 
surements differ by up to a few tenths from the Tycho magnitudes measured in 
the early Nineties. We could measure the short-term (rotational) variabilities 
of TYC 9300-0891-lAB and of TYC 9300-0529-1, which has amplitude up to 
~0.10 mag in the V bancp~| (see Sect. 4). For the long-term variability, we only 
have integrated magnitudes (owing to the low angular resolution) of all three 


components provided by the All Sky Automated Survey (ASAS; Pojmanski 


1997) (see Fig.[^, which show the effect resulting from the combination of 


the brightness variability of the three stars. Since all components have similar 


Inferred from R-band light curve amplitudes (Messina et al. 2014) 
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brightness, in principle they all can signihcantly contribute to the observed 
variability. 

The Li abundance, X-ray emission, photometric variability, and fast rotation 
are all consistent with the inferred young age of this system. 


3 Observations 


This stellar system was observed in 2011 at the South African Astronomi¬ 
cal Observatory (SAAO). We could observe it for 6 nights, from April 3 to 
12 (see Table[^. We used the 1-m telescope equipped with STE4 CCD cam¬ 
era and Johnson-Cousins BV(RI)c hlter set. We could collect a total of 31 
frames in the R band using integration time of 10 sec. The data reduction 
was performed using the IRAFp] tasks within DAOPHOT. All frames were 
bias subtracted and flat held corrected. We used aperture photometry to ex¬ 
tract the series of magnitudes of all stars detected in the 5'x5' held of view. 
The two stars 2MASS J18501087-7156309 (C) and 2MASS J18495139-7155131 
(CK) (see Table§, served as comparison and check, respectively, to get the 
diherential magnitude time series. The achieved photometric precisions of the 
diherential magnitudes were 0.002 mag and 0.003 mag for TYC 9300-0891-1 
and TYC 9300-0529-1, respectively. 

We hnd the magnitude of the comparison to stay constant during the run with 
a standard deviation ac-CK = 0.013 mag. The two components A and B of 
TYC 9300-0891-1 in almost all nights could not be spatially resolved, and we 
extracted their integrated magnitudes series. 


The SAAO STE4 camera mounted on the Im telescope has a plate scale of 
0.31"/px. On few best-seeing nights, we could resolve partially component A 
from B, but sufficiently to extract accurate magnitudes with the PSF tech¬ 
nique. In the left panel of Fig.[^ we show a section of a frame centred on 
the stellar system and collected on April 10, 2011 where the binary nature of 
TYC 9300-0891-1 is clearly visible. In the right panel, we plot the surface plot 
obtained with IRAF which shows the A and B component partially resolved. 
The PSF photometry allowed us to extract the following magnitudes Va= 
11.65±0.011, Vb= 11.67±0.010, and Vc= 11.68T0.008 mag. We measured a 
separation p = 0.95±0.47" between the A and B components, in good agree¬ 
ment with the separation measured by Tycho, and similar position angle PA 
= 130°±10 (or PA = 310°±10 depending on which component either A or B 

^ IRAF is distributed by the National Optical Astronomy Observatory, which is 
operated by the Association of the Universities for Research in Astronomy, inc. 
(AURA) under cooperative agreement with the National Science Foundation. 
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is assumed to be the primary). 


We carried out a second observing run in 2013 for 45 nights from May 31 until 
October 9 (see Tabl^ at the Klein Karoo Observatory (225 m a.s.h, Western 
Cape, South Africa) with a 30-cm (f/8) RCX-400 telescope having 21'xl4' 
held of view. It is equipped with SBIG ST8-XME CCD camera and Johnson 
V and Schuler R hlters. We could collect a total of 1034 frames in the R band 
using integration time of 27 sec. In Nov 21, 2014 we could collect 7 additional 
frames in the V hlter. On each night, from 4 to 10 consecutive frames were 
summed up to improve the S/N ratio to get on average from one to three mag¬ 
nitudes per night to be used for the subsequent analysis. The data reduction 
was performed as for the hrst run data set collected in 2011. Despite the smaller 
telescope aperture, TYC 9300-0891-1 and TYC 9300-0529-1 were spatially well 
resolved by our observing system, then we used again aperture photometry to 
extract the series of magnitudes. The best extraction aperture was changed on 
each night according to variable seeing conditions. Owing to the larger FoV, 
we could identify a new set of comparison stars (see Tablej^ with more sta¬ 
ble light curves. The star TYC 9300-0573-1 (C) was used as comparison star 
to get differential values of our targets. The achieved photometric precisions 
of the differential magnitudes were 0.004 mag and 0.005 mag for TYC 9300- 
0891-1 and TYC 9300-0529-1, respectively. The two stars TYC 9300-1261-1 
(CKl) and 2MASS J18485976-7152445 (CK2) served as check stars. We hnd 
the magnitude of the comparison to stay constant during the whole run with 
standard deviation <Jc-cki = 0.009 and <Jc-ck2 = 0.011 mag. Of course, the 
two components A and B of TYC 9300-0891-1 could not be spatially resolved, 
then we extracted their integrated magnitudes series. 


4 Periodogram analysis 


Our targets, being young low-mass stars, are expected to host some level of 
magnetic activity and, therefore, are expected to exhibit starspots on their 
photospheres. The light modulation induced by the star’s rotation, which car¬ 
ries these spots in and out of view, is expected to be quasi-periodic and to have 
the star’s rotation period. In these spotted stars, the rotation period can be 
derived by Fourier analysis of the magnitude time series. In our specihc case, to 


search for the photometric rotation periods we used Lomb-Scargle (LS; Scargle 


1982) and Clean (Roberts et al.j~987|) periodograms. An estimate of the False 


Alarm Probability (FAR), that is the probability that a peak of given power 
in the periodogram is caused by statistical variations, i.e., by Gaussian noise, 
was done using Monte Carlo simulations according to the approach outlined 


by Herbst et ah (2002). A more detailed description is given by Messina et ah 


(2010). The uncertainty on the rotation period determination was estimated 
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following Lamm et al. (2004) (see also Messina et al. 2010). 


In the case of TYC 9300-0891-lAB, as already mentioned, the components A 
and B are unr esolved. However, since they have same brightness, in principle 
both components may signihcantly contribute to the observed variability leav¬ 
ing signs of their periodicities (that is two power peaks) in the periodogram. 
The LS periodogram in both seasons revealed three signihcant periodicities 
(FAP < 1% in both seasons): the dominant period is P = 1.383±0.001 d, the 
secondary periods are P = 3.60±0.01 d, and P = 0.560 d (see top of Fig.[^. The 
Clean periodogram showed P = 1.383 d, whereas the longest and the shortest 
periods were effectively hltered out by the Clean algorithm. P = 1.383 d is 
the stellar rotation period of either component A or B (more specihcally, of 
the component exhibiting the larger variability), whereas P = 3.60 d and P = 
0.560 d are beat periods between P = 1.383 d and the data sampling of about 
Id, imposed by the rotation of the Earth and the hxed longitude of the obser¬ 
vation site. The light curve amplitude (amplitude of the htting sinusoid) is AR 
= 0.02 mag. The presence of only one power peak in the cleaned periodogram 
indicates that only one component is effectively contributing to the observed 
variability, the other component having either no spots at the epochs of our 
observations or a quite uniform spot distribution in longitude that produces 
very small rotational modulation, or spots out of view (see Sect. 5.1). In the 
top panels of Fig.[^ we plot the KKO differential magnitude timeseries; the 
LS periodogram, where the red dotted line is the spectral window and the 
horizontal dashed line the power corresponding to 99% conhdence level; the 
Clean periodogram; and the light curve phased with the rotation period with 
overplotted a sinusoidal ht (diamonds are SAAO data, dots are KKO data). 


In the case of TYC 9300-0529-1, our results are summarized in Fig.|^ The 
LS periodogram in both seasons reveals again three signihcant periodicities: 
the dominant period is P = 1.634±0.001 d, and the secondary periods are P 
= 2.57±0.01d, which is the same period reported by Messina et al. (2011) 
and derived from the analysis of ASAS timeseries, and P = 0.62 d. The Clean 
periodogram reveals the longer periods in the hrst run, but only the shorter 
period P = 1.632 d in the second run. All secondary periods are again alias 
arising from the spectral window and are effectively hltered out by the Clean 
algorithm. We conclude that P = 1.634 d is the stellar rotation period. The 
observed light curve has amplitude AR = 0.08 mag in the hrst season. In the 
second season, we note that the amplitude has progressively increased from 
AR = 0.03 mag up to AR = 0.06 mag, as as consequence of the active regions 
evolution. 


In the left panel of Fig.|^ we plot the distribution of rotation periods versus 
B—V colors of Octans members, which is updated with respect to [Messina et 












al. (2011). We find that our targets have their rotation periods well within 


the distribution of the rotation periods of other Octans members in the same 
color range (0.6 < B—V < 0.8). This circumstance further supports the mem¬ 
bership of our target to the Association. In the right panel of Fig.[^ we plot 
the distribution of V-band light curve amplitudes of Octans members versus 
rotation periods. In our case, the R-band light curve amplitudes of our targets 
were multiplied by a factor 1.2 to be comparable to the V-band amplitudes 
( Messina et al.||2014 ). TYC 9300-0891-1 is positioned in the lower boundary of 
the amplitude distribution, owing to the very low inclination of the rotation 
axis (see Sect. 5.1). 


5 Target’s parameters 

5.1 TYC 9300-0891-1A and B 


We can use effective temperature and luminosity to derive the stellar radius 
and the inclination of the rotation axis. Assuming d = 173 pc (see Sect. 2), 
a mean value of the observed magnitude V = 11.66 mag for both compo¬ 
nents, bolometric correction BC = —0.33mag ( |Pecaut fc Mamajek 2013), 
Teff = 5168K (da Silva et al. 2009), we infer the bolometric magnitude Mboi 
= 5.13±0.05mag, luminosity L = 0.70±0.06Lq, and the stellar radius R = 
1.O4TO.12R0. Combining the radius with the projected rotational velocity 
and the rotation period P = 1.383 d, we infer the inclination of the stellar 
rotation axis i ~ 14° for either component A or B. 

To check whether TYC 9300-0891-1 could generate detectable flux rotational 
modulation despite such very inclined rotation axis, we modeled the observed 
KKO R-band light curve using Binary Maker V2.0 (Bradstreet 1993). Bi¬ 
nary Maker V2.0 models are almost identical to those generated by Wilson- 
Devinney program (Wilson & Devinney 1971) and uses Roche equipotentials 
to create star surfaces. In our modeling the second component is essentially 
“turned off” (i.e. assigned a near zero mass and luminosity) in order to model 
a single rotating star (see Messina et al.||2015 ). In Sect. 4, we showed evidence 
that only one component has produced the observed variability. One circular, 
cool spot of uniform temperature was assumed and the light curve was £t by 
manual iteration. Gravity-darkening coefficient has been assumed u = 0.25 
(Kopal 1959), and limb-darkening coefficient for the R-band {ur = 0.60) of 


Al-Naimy (1978) was also adopted. We adopted the effective temperature of 
the star Teg and inclination i = 14° as input parameters. The temperature 
contrast between spots and surrounding photosphere was hxed to T spot /T phot 
= 0.8, which is a value generally measured in KOV spotted stars (see, e.g.. 


Berdyugina 2005). The ’dilution’ effect in the amplitude of the light curve 


(in our case by a factor 2) arising from the constant flux contribution from 
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the unresolved companion component was also considered (using the ’third 
light’ option within the model input parameters). The only free parameters 
in our model remained the spot hlling factor and the spot latitude/longitude. 
Our aim is not to hnd unique solution of the spot area and position, but to 
see if conhgurations exist that can allow rotational modulation of 0.02 mag 
amplitude. In Fig.[^ we plot as an example one of such solutions where in the 
top we plot the light curve (crosses represent the normalized flux values) with 
overplotted the model £t, and in the bottom we plot the spot configurations 
corresponding to the minimum and maximum spot visibility at phases 0 = 
0.19 and 0 = 0.69, respectively. In the specihc example, we used one circular 
spot with 25° radius, positioned at the equator. As we can see, despite the 
very low inclination value, the star can modulate the spots visibility giving 
rise to low amplitude flux rotational modulation. 

It is likely that at the epochs of our photometric monitoring, one of the two 
components (either A or B) had either no spots or spots, but out of view owing 
to the low inclination, thus giving no contribution to the observed variability 
and preventing us from the rotation period determination. 

We use the available optical, near-IR, and IR photometry to build the observed 


Spectral Energy Distribution (SED). The BVI magnitudes from Torres et al. 


(2006), JHK magnitudes from 2MASS ( Cutri et'ah]|2003[ ), and W1-W4 mag¬ 
nitudes from WISE (Cutri et al. 2003) are listed in Table|^ We note that all 
magnitudes refer to the unresolved AB components. The SED was htted with 
a grid of theoretical spectra from the NextGen Model ( [Allard et al.||2012 ) and 
the best £t is obtained with a model of Te^ = 5100T100 K, log g = 4.0±0.3 dex 
and metallicity [Fe/H] = 0.0 (see top panel of Figj^. We note no evidence of 
IR excess. 


5.2 TYC 9300-0529-1 


Assuming d = 175 pc (see Sect. 2), brightest observed magnitude V = 11.68 
mag, bolometric correction BC = —0.33mag (Pecaut & Mamajek 2013) , Teff 
= 5140K (da Silva et al. 2009), we infer the bolometric magnitude Mboi = 
5.13±0.05 mag and luminosity L = O.69±O.O6L0, and the stellar radius R = 
1.05±0.11 Rq. From these values we infer the inclination of the stellar rota¬ 
tion axis i ~ 50°. 

Although TYC 9300-0529-1 and TYC 9300-0891-1 have same masses, simi¬ 
lar rotation rates, and similar activity levels, as expected from the empirical 
rot at ion-activity relations ( Messina et al.||2011 ), the larger inclination of TYC 
9300-0529-1 rotation axis allows a larger amplitude rotational modulation (see 
Fig.|4|). 

Also for this component we use the available photometry to build the observed 
SED (see also Table[^. The SED was htted with a grid of theoretical spectra 
from the NextGen Model (Allard et al.||2012) and the best £t is obtained with 
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a model of Teff = 5100±100 K, log g = 4.0±0.3 dex and metallicity [Fe/H] = 
0.0 (see bottom panel of Figj^. Also for this star we note no evidence of IR 
excess. 


To infer the mass and pnt some constrain on the age of the system’s compo¬ 
nents, we compared Teg and Inminosity with a set of evolntionary tracks and 
isochronones taken from Baraffe et al. (1998) for solar metallicity. In Fig.|^ 
we see that all three components are £t by an evolntionary track of M = 
0.95±0.05Mq, and by an isochrone of 35±5 Myr, if we take into acconnt the 
uncertainties. This result further supports the coevalness of all components. 


6 Discussion 


At an age of about 35 Myr, our targets should have passed the disc-locking 
phase, and should be experiencing the phase of stellar radius contraction (see, 
e.g., Bouvier et al. 2014). Owing to angular momentum conservation, their 


rotation rates are increasing while approaching the Zero-Age Main Sequence. 

Differently than the mentioned case of BD—211074, whose components have 
rotation periods differing by about 70%, in the stellar system under analysis 
the difference amount to only 16%. Considering that masses, ages, and ini¬ 
tial chemical compositions are almost identical for all components, only disc 
lifetimes and initial rotation periods may cause this small difference in the 
rotational evolution of the components. TYC 9300-0891-lAB is a close binary 
where the component A has a companion B at a distance of about 157 AU. We 
want to probe whether this difference in rotation periods may be due to the 
dynamical effects of component B on the disc of component A (or viceversa) 
that have facilitated its dispersal, making shorter the disc lifetime and, con¬ 
sequently, the disc-locking phase. As a consequence, TYC 9300-0891-1 should 
have started experiencing the rotation spin up earlier than the more distant 
component TYC 9300-0529-1, determining the observed rotation period dif¬ 
ference. 

The accretion disc response to an external perturber has been already subject 


of studies. Clarke &: Pringle (1993) analysed the effect of disc dispersal in the 


case of stellar fly-by. Depending on the kind of encounter, coplanar prograde, 
coplanar retrograde, and orthogonal, the effect will be more or less destructive 
on the disc life, the hrst conhguration being the most destructive. Similar re¬ 


sults are found by Heller (1995) who computes that less than about half of the 


initial disc mass is lost during one single encounter, which may be captured 
by the perturber or ejected by the system. A case more similar to the one 


under our analysis is discussed by Korycansky et al. (1995). They show that 


it is possible to model the accretion disc response in an eccentric binary sys- 
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tern with succession of parabolic passages. Then, the disc undergoes a series 
of impacts and looses angular momentum to the secondary, thus reducing its 
lifetime. 

To check our hypothesis, we use as our guideline the angular momentum evolu¬ 
tion model of Gallet & Bouvier (2013). We see in their Fig. 3 that our targets 
are mostly similar to the fast rotators modeled by them. In the age range 
from about 5 to 30 Myr, the angular velocity linearly increases in log(G/ff 0 )- 
log(age) scale, allowing us to estimate the disc lifetime difference corresponding 
to the rotation period difference between TYC 9300-0891-1 (hereafter stari) 
and TYC 9300-0529-1 (hereafter star 2 ). 

However, we hrst need to make some assumption on the initial values of the 
rotation periods of stari and star 2 during the disc-locking phase. 

The rotation period distribution of the Orion Nebula Cluster members at 
about 1 Myr can be considered as representative of the period distribution 
during the disc-locking phase. We thus retrieved from Herbst et al.| (2002) 
the rotation periods within a small mass range, 0.8-1.2 Mq, bracketing our 
components. The period distribution is found to be double peaked with the 
shorter periods ranging from about 1 to 5 days. Our hrst assumption is that 
our stari and star 2 have their within this range during the disc locking 
phase. A second assumption is that after 10 Myr the disc-locking in no longer 
effective and stars are free to spin up. 

Based on these assumptions, on the present rotation periods, and on the rate 
of spin up in the 10-30 Myr interval from the Callet & Bouvier model, we 
derive that stari during the disc-locking phase can have initial rotation period 
in the range from 5 to 2 days and star 2 in the range from 5 to 3 days (see left 
panel of Fig.[^. Comparing these two relations we can compute the possible 
combinations of initial values that will produce the presently observed rotation 
periods (see rectangular area in the right panel of Fig.|^. More specihcally, 
the solid line represents the combinations when we assume both stars to have 
same disc lifetimes, whereas the upper region corresponds to the assumption 
that Ti < T 2 , and the lower region to the assumption ri > r 2 , where ri and 
T 2 are the disc lifetimes of stari and star 2 . The case ri < T 2 has a probability 
of 57% with respect to the opposite case. However, on the light of the uncer¬ 
tainties in the comparison with the model, such a difference in probability is 
not signihcant and our hypothesis of a lifetime shortening of stari by dynami¬ 
cal effect by the perturber, although possible, is not favored (more probable). 
Both scenarios t\ < T 2 and ri > T 2 can account for the observed rotation 
period difference. 

A plausible scenario that comes out from our investigation is that the system’s 
components sheared similar initial rotation periods, and during the 35 Myr of 
life have undergone similar rotation period evolution being about similar their 
disc lifetimes. Therefore, in this specihc system it seems that a stellar com¬ 
panion at a distance of about 157 AU may be still distant enough to produce 
no effective perturbation on the disc of the other star. The perturbation, if 
present, but our analysis does not secure us on this issue, may only slightly 
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shorten the disc lifetime to cause a 16% difference, at most, between the ro¬ 
tation periods. 


On the contrary, in the mentioned visual triple system BD—211074 (Messina 


et al.|20i4) in the 20-Myr (3 Pictoris association, where the component BD—21 


1074B has a perturber companion BD—211074C at only 15 AU of distance, 
whereas the component BD—211074A is at 160 AU, the enhanced disc dis¬ 
persal was effective and allowed the secondary component B to start its spin- 
up at much earlier stages, reaching a rotation period about 70% faster than 
BD-211074A. 


7 Conclusions 


We have collected from the literature information that allows us to conh- 
dently state that the components TYC 9300-0891-lAB -|- TYC 9300-0529-1 
are physically bound forming a triple system and they are member of the Oc- 
tans association. We have carried out a photometric monitoring of this system 
in 2011 and 2013. We could determine the rotation period P = 1.383d for 
the unresolved components TYC 9300-0891-lAB and P = 1.641d for TYC 
9300-0529-1. Using archival data, we have modeled their SEDs and we did 
not hnd any evidence for IR excess. We derive the luminosity of all compo¬ 
nents and comparing with evolutionary model we inferred for the system an 
age of 35±5 Myr, in good agreement with their membership with the Octans 
Association. Modeling the photometric variability, we found that the photo- 
spheric activity level is signihcant with cool spots covering at least 10% of 
the stellar disc. Since all three components have almost identical mass, age, 
and initial chemical composition, the rotation period difference likely arises 
from either different values of the initial rotation periods or different lifetimes 
of their discs. However, the 16% rotation period difference between the com¬ 
ponents is too small to infer meaningful information on possible differences 
in the disc evolution history by a comparison with the angular momentum 
evolution model of Gallet & Bouvier (2013). 

In the previously examined case of BD—211074 in the f3 Pictoris Association, 
the period difference was sufficiently large (70%) to infer that the closer com¬ 
panion at 15 AU determined a more efficient disc dispersal. 

We are going to study other similar systems to see on a statistical basis how 
much signihcant is the role played by such perturber in the primordial disc 
lifetime or if the observed period distribution is dominated by different values 
of the initial rotation periods. 
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Fig. 1. V-band magnitude timeseries of the spatially unresolved system TYC9300 
0891 lAB + TYC9300 0529 1 as observed by ASAS and at KKO in the 2014 (open 
square). 
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Fig. 2. top panel: section of a frame collected with the Im SAAO telescope. 
The components A and B are at bottom left, bottom panel: Surface plot of 
TYC 9300-0891-1AB where the two components are partially resolved. 
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TYC9300-0891-1A P3 


1 .381 


Lomb —Scorgle (KKO data) 


Cleon (KKO dato) 



Fig. 3. Top panels: Results of periodogram analysis of TYC 9300-0891-lA. We 
plot the R-band magnitude timeseries versus Heliocentric Julian Day; the Lomb-S- 
cargle periodogram (solid line) with overplotted the spectral window (dotted red 
line), whereas the horizontal dashed line represent the 1% FAR level; the Clean 
periodogram. Bottom panel: the phase light curve with the rotation period and 
overplotted the sinusoidal fit (solid red line). 
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Fig. 4. Results of periodogram analysis TYC9300-0529-1 (see Fig.|^for explanation). 
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Fig. 5. Distribution of rotation periods versus B—V color (left) and light curve 
amplitudes versus period (right) of Octans members. 
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Fig. 6. Results of spot modeling of TYC 9300-0891-1. top: the normalized light curve 
(crosses) with amplitude AF/F = 0.037 phased with the P = 1.383 d rotation period 
and overplotted the model fit (solid line), bottom: spot configuration at minimum 
and maximum spot visibility. 
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10“ lO'^ 

A (A) 


Fig. 7. Spectral Energy Distributions of TYC 9300-0891-lAB (top) and 
TYC 9300-0526-1 (bottom). Red bullets are the observed fluxes (see correspond¬ 
ing magnitudes in Tablej^ whereas blue lines are the best fit synthetic spectra from 
NexGen Models. The error associated to the fluxes are smaller than the symbol size. 



Fig. 8. HR diagram for TYC 9300-0891-1 and TYC 9300-0529-1. The dashed line is 
the isochrone corresponding to an age of 35 Myr and the solid line is the evolutionary 
mass track corresponding to 0.95 Mq from Baraffe et al. (1998). 
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Initial rotation period (d) 


Initial period star 2 (d) 


Fig. 9. Left panel: possible values of disc lifetimes versus initial rotation period 
according to the Gallet &: Bouvier (2013) model for TYC 9300-0891-1 (stari) and 
TYC 9300-0529-1 (star 2 ). Right panel: the rectangular area limits the possible com¬ 
binations of initial rotation periods according to the mentioned model. The currently 
observed rotation periods can be obtained assuming either ti < T 2 , or ri > T 2 , or 
n = t 2 . 
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Table 1 

Observation log with listed dates and number of frames in the R filter. 


2011/04/03 

3 

SAAO 

2011/04/05 

4 

2011/04/09 

6 

2011/04/04 

5 

2011/04/08 

5 

2011/04/12 

8 

2013/05/31 

8 

KKO 

2013/07/12 

6 

2013/08/13 

22 

2013/06/05 

12 

2013/07/15 

5 

2013/08/23 

14 

2013/06/09 

13 

2013/07/17 

14 

2013/08/24 

9 

2013/06/12 

13 

2013/07/21 

6 

2013/08/25 

6 

2013/06/15 

243 

2013/07/22 

12 

2013/08/26 

17 

2013/06/16 

104 

2013/07/24 

13 

2013/09/01 

21 

2013/06/17 

53 

2013/07/26 

5 

2013/09/02 

14 

2013/06/19 

18 

2013/07/28 

11 

2013/09/03 

6 

2013/06/20 

10 

2013/07/29 

7 

2013/09/04 

17 

2013/06/26 

84 

2013/07/30 

13 

2013/09/05 

14 

2013/07/06 

19 

2013/07/31 

24 

2013/09/06 

8 

2013/07/07 

7 

2013/08/01 

28 

2013/09/12 

6 

2013/07/08 

6 

2013/08/02 

19 

2013/09/26 

16 

2013/07/09 

8 

2013/08/03 

6 

2013/09/29 

7 

2013/07/10 

17 

2013/08/05 

21 

2013/10/07 

6 

2013/07/11 

23 

2013/08/10 

15 

2013/10/09 

8 





2014/11/21 

7(V) 
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Table 2 

Comparison stars 


name 

RA(J2000) DEC(J2000) 

Mag 


SAAO 


2MASS J18501087-7156309 (C) 

18 50 10.88 -71 56 30.9 

J = 12.16 

2MASS J18495139-7155131 (CK) 

18 49 51.40 -71 55 13.2 

J = 11.82 


KKO 


TYC 9300-0573-1 (C) 

18 48 29.114 -71 53 22.48 

V = 11.6 

TYC 9300-1261-1 (CKl) 

18 48 23.246 -71 57 15.55 

V = 11.98 

2MASS J18485976-7152445 (CK2) 

18 48 59.763 -71 52 44.56 

J = 11.18 


Table 3 

Magnitudes from the literature. 



B 

V 

I 

J 

H 

K 

W1 

W2 

W3 

W4 

TYC 9300-0891-lAB 

11.63 

10.86 

9.96 

9.31 

8.89 

8.73 

8.65 

8.67 

8.60 

8.29 

TYC 9300-0529-1 

12.39 

11.59 

10.68 

10.04 

9.63 

9.52 

9.48 

9.48 

9.45 

<8.93 
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Table 4 

Summary of information on TYC 9300-0891-lAB and TYC 9300-0529-1 from the 
literature and the present study. 


TYC9300-0891-1AB 

TYC9300-0891-1A TYC9300-0891-1B TYC 9300-0529-1 


Sp. Type 


KOVe 

KOV 

KOV 


V (mag) 


11.65 

11.67 

11.68 

this study 

Mass (Mq) 


0.95±0.05 

0.95±0.05 

0.95±0.05M 

this study 

Teff (K) 


5168±100 


5140±100 

D09 

Bum (Lq) 


0.70±0.06 

0.70±0.06 

0.69±0.06 

this study 

Radius (R©) 


1.04±0.12 

1.04±0.12 

1.05±0.12 

this study 

log g 


4.0±0.3 


4.0±0.3 

this study 

[Ee/H] 


0.0 


0.0 

this study 

EW Li (mA) 

310 


300 

T06 

age (Myr) 


35±5 

35±5 

35±5 

this study 

UVW (kms- 


-13.7; -3.5; - 

11.5 

-13.3;-4.4;-11.4 

M14 

XYZ (pc) 


124.7; -93.9; - 

-74.7 

126.1; -95.0; -75.5 

It? ?? 

masyr“^ 

(Tycho) 

14.8±2.1 


14.9±2.1 

ESA (1997) 

fis masyr“^ 

(Tycho) 

-16.1±2.0 


-12.3±1.9 


fia masyr“^ 

(WDS) 

15 


14 

MOl 

fj,s masyr“^ 

(WDS) 

-16 


-18 


fj,a masyr“^ 

(UCAC2) 

14.4±5.8 


14.3±1.5 

Z04 

fis masyr“^ 

(UCAC2) 

-20.4±5.8 


-18.4±3. 


masyr“^ 

(PPMXL) 

14.8±2.1 


14.1±1.6 

RIO 

US masyr“^ 

(PPMXL 

-16.1±2.0 


-16.4±1.6 


V,. (kms“^) 


-3.0 


-2.1 

T06 





-2.26±0.47 

E14 

distance (pc 

) 

173 


175 

M14 


TOO: Torres et al. (2006); M14: Murphy et al. (2014); Z04: Zacharias et al. (2004); 
RIO: Roeser et al. (2010); E14: Elliott et al. (2014); MOl: Mason et al. (2001); 
D09: da Silva et al. (2009). 
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